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Gamma-ray burst (GRB) dust echoes were first proposed as an alternative 
explanation for the supernova-like (SN-like) components to the afterglows of GRB 
980326 and GRB 970228. However, the spectroscopic identification of Type Ic SN 
2003dh associated with GRB 030329, as well as the identification of SN-like components 
to the afterglows of other GRBs, appears to have confirmed the GRB/SN paradigm. 
However, the likely progenitors of Type Ic SNe are Wolf-Rayet WC stars, and late-type 
WC stars have been observed to be surrounded by dust, at a distance of 10 14 - 10 15 cm 
■ from the star. Consequently, we revisit the possibility of GRB dust echoes, not on a 

• timescale of weeks after the burst but on a timescale of minutes to hours. We find that 

if the optical flash is sufficiently bright and the jet sufficiently wide, GRB afterglows 
may be accompanied by chromatic variations on this timescale. From these signatures, 
such model parameters as the inner radius of the dust distribution, the initial opening 
Q-f angle of the jet, etc., may be deduced. With rapid and regular localizations of GRBs 

by HETE-2, Integral, and now Swift, and new and improved robotic telescope systems, 
these early-time GRB dust echoes may soon be detected. We describe one such robotic 
telescope system, called PROMPT, that the University of North Carolina at Chapel 



Hill is building at the Cerro Tololo Inter- American Observatory in greater detail. 



Subject headings: dust, extinction — gamma rays: bursts — stars: mass loss — stars: 
Wolf-Rayet — supernovae: general — telescopes 



1. Introduction 

GRB dust echoes were first proposed by Waxman & Draine (2000) and Esin & Blandford 
(2000) as an alternative explanation for the SN-like components to the afterglows of GRB 980326 
(Bloom et al. 1999) and GRB 970228 (Reichart 1999; Galama et al 2000). Here, "dust echo" 
refers to light that is either (1) scattered by dust into the line of sight or (2) absorbed by dust 
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and then thermally reemitted into the line of sight. In both cases, light is delayed due to the 
greater path lengths. Reichart (2001) modeled and computed dust echo light curves and spectral 
flux distributions (SFDs) for both cases and found that while dust echoes can mimic SN light 
curves they cannot mimic SN SFDs, at least not near the spectral peak: For example, the second 
component to the afterglow of GRB 970228 cannot be explained by a dust echo because its 
spectral peak is too narrow for the first case and at too high of a frequency for the second case. 

The spectroscopic identification of SN 2003dh associated with GRB 030329 (Stanek et al. 
2003), as well as the spectroscopic identification of SN 1998bw associated with the unusual GRB 
980425 (Galama et al. 1998) and the identification of SN-like components to the afterglows of 
at least seven other GRBs (e.g., Bloom et al. 1999; Reichart 1999; Galama et al. 2000; Bloom 
et al. 2002; Garnavich et al. 2003; Price et al. 2003; Delia Valle et al. 2003; Cobb et al. 
2004; Levan et al. 2005), appears to have confirmed the GRB/SN paradigm, at least for most 
long-duration/soft-spectrum GRBs. 

The leading models for making GRBs from collapsing massive stars are the collapsar model 
(e.g., Woosley 1993) and the supranova model (e.g., Vietri & Stella 1999), both of which result in 
the formation of a black hole and an accretion disk from which an ultra-relativistic jet is produced, 
but in the case of the supranova model a neutron star that must first shed angular momentum 
forms first. In the case of the collapsar model, the progenitor must first shed its outer layers 
in a Wolf-Rayet (WR) wind (or be stripped of them by a companion) if the jet is to escape the 
progenitor. Indeed, both SN 2003dh and SN 1998bw are Type Ic SNe, the progenitors of which 
have lost both their hydrogen and helium layers. 

Given that the SN start time matches the GRB time to within about two days for both of 
these events (e.g., Iwamoto et al. 1998; Hjorth et al. 2003), the collapsar model is favored, at 
least for these events. However, given that there is no overlap between the sample of GRBs for 
which SNe and SN-like components to afterglows have been detected and the sample of GRBs for 
which blueshifted soft X-ray lines have been detected (e.g., Piro et al. 2000; Reeves et al. 2002, 
2003; Butler et al. 2003; Watson et al. 2003) - which can be interpreted in terms of the supranova 
model - a dichotomous scenario cannot yet be ruled out. However, in this paper we will narrow 
our focus and consider only GRBs with WC (or WO) progenitors, the most evolved of the WR 
stars. 

Given that late- type WC stars are expected to be surrounded by dust at a distance of >10 13 
- 10 14 cm from the star and have been observed to be surrounded by dust at a distance of 10 14 - 
10 15 cm - not at a distance ~10 18 cm as would be required to mimic a SN light curve - we now 
revisit the possibility of GRB dust echoes, not on a timescale of weeks after the burst but on a 
timescale of minutes to hours. In §2, we summarize both the expectation for and observations of 
dust at these distances from late-type WC stars. In §3, we present model dust echo light curves 
and color histories for a range of parameter values and discuss limitations of the model. In §4, we 
draw conclusions and summarize a robotic telescope system that should be ideal for identifying 
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and studying dust echoes at early times after the burst. 



2. Dust around Late- Type WC Stars 

Many studies of infrared emission from late-type WC stars have shown them to be dust-making 
machines (e.g., Allen, Harvey & Swings 1972; Gehrz & Hackwell 1974; Cohen, Barlow & Kuhi 
1975; Williams, van der Hucht & The 1987). However, grains cannot condense out of the dense 
stellar wind within the star's dust sublimation radius r s . For an optically thin medium in radiative 
equilibrium: 

4vra 2 Q(^T>T 4 = ?ra 2 Q(^^^, (1) 



where na 2 Q(a, T) is the mean absorption cross-section for a grain of size a averaged over a thermal 
spectrum of temperature T, T s is the sublimation temperature, and T* and R* are the effective 
temperature and radius of the star. Taking Q(a,T s ) ~ Q(a, T*) (e.g., Tuthill, Monnier & Danchi 
2001), simplification yields: 

1 tT,^ 



2 KT.J R " (2) 

For an effective temperature T* ~ 35,000 K and radius ~ 2 x 10 11 cm (e.g., Crowther 1997), and 
a sublimation temperature T s ~ 1,500 K (e.g., Tuthill, Monnier & Danchi 2001), the sublimation 
radius is r s ~ 5 x 10 13 cm. This is of course only a rough estimate: There is much ambiguity in 
how one measures the effective temperature and radius of a WR star because of the opacity of its 
wind, and the sublimation temperature varies with grain composition. Furthermore, the process 
by which dust forms beyond the sublimation radius is not fully understood: Dust is not expected 
to be present in large quantities at this innermost possible radius. 

Observations of dust around late-type WC stars yield distances that are consistent with 
this expectation. Williams, van der Hucht & The (1987) modeled infrared SFDs of 23 WC stars 
assuming an isotropic distribution of dust and found inner radii that are typically in the range 
10 14 - 10 15 cm. Others have measured the dust distribution directly using one-dimensional speckle 
interferometry (Allen, Barton & Wallace 1981; Dyck, Simon & Wolstencroft 1984 for WC9 star 
WR 104), aperture masking with Keck (Tuthill, Monnier & Danchi 1999, 2002 for WR 104), lunar 
occultation (Ragland & Richini 1999 for WC9 star WR 112; Mondal & Chandrasekhar 2002 for 
WR 104); and direct imaging with HST-NICMOS2 (Marchenko, Moffat & Grosdidier 1999 for 
WC7 star WR 137). However, all of these might be special cases: WR 104 and WR 112 are 
the two most extreme dust producers in the sample of Williams, van der Hucht & The (1987), 
presumably due to wind-wind collisions with short-period companions, which additionally result 
in a highly anisotropic dust distribution; and WR 137, which is technically not a late-type WC 
star, is an episodic dust producer due to wind-wind collisions with a long-period companion on an 
elliptical orbit. Despite these caveats, typical distances to their dust distributions are 4 x 10 14 cm 
for WR 104, 6 x 10 14 cm for WR 112, and 7 x 10 15 cm for the less relevant WR 137. 
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3. GRB Dust Echo Model and Limitations 

We now present dust echo light curves and color histories using the model of Reichart (2001), 
but for significantly smaller dust distributions than he explored. We summarize the dust and light 
geometries and model parameters in the top panel of Figure 1: Dust does not form within radius 
R of the progenitor and is destroyed within half-angle 9j e t of the jet axis by bipolar jets (e.g., 
Waxman &i Draine 2000; Fruchter, Krolik & Rhoads 2001). Exterior to the solid curve, the dust 
density n(r) oc r~ 2 . The optical depth through this dust distribution is t v = n{R)a u R, where a v 
is the total (absorption plus scattering) cross section of the dust grains at frequency v. A pulse of 
light is assumed to be emitted within half-angle 9 co \ of the jet axis (Reichart 2001). 

Light that is emitted at angle 9 relative to the line of sight, where 9j et < 9 < 6 co i, and 
is scattered by dust at radius r > R into the line of sight will arrive at the observer with an 
observer-frame time delay given by: 

t= r(1 + ^ (l-cos 0). (3) 

In the bottom panel of Figure 1 we plot curves of constant arrival time for an on-axis jet. The 
spectral flux of the dust echo at time t is then given by integrating along the illuminated portions 
of the corresponding curve the product of (1) the spectral flux of the pulse of light, (2) the dust 
density n(r,6), (3) the probability that the light will scatter an angle 9 into the line of sight, 
and (4) e~ Tl, ( 1+z ^ r,e \ where t„(i +z ) (r, 0) is the optical depth integrated along the light path. 
Expressions for the differential scattering cross section and TW 1+2 .)(r, 0), as well as useful analytic 
approximations for the light curve and color history of the dust echo, can be found in §2.1 of 
Reichart (2001). 

In Figures 2 - 5, we plot dust echo light curves and color histories in which we vary each of 
the four model parameters - R, 9j e t, t^( 1+2 ), and 9 co i - while holding the other three constant: 

Figure 2: Greater values of R result in later turn-on times because of the greater path length that 
the light must travel. Greater values of z result in later turn-on times because of cosmological 
time dilation. The dust echo starts out bluer than the optical flash because blue light scatters 
preferentially, but grows redder than the optical flash because of increasing path lengths through 
dust. The re-brightening of the dust echo at later times corresponds to light from the jet that is 
pointed away from us backscattering off of the far side of the dust distribution. This light is bluer 
because destruction of dust by the two jets leaves little dust in the path of this light to redden it. 

Figure 3: Greater values of 9j e t result in later turn-on times because of the greater path length 
that the light must travel. The re-brightening of the dust echo at later times occurs an equal 
amount of time earlier. 

Figure 4-' Greater values of T v n +Z \ result in dust echoes that grow redder because of increasing 
path lengths through greater dust densities. 

Figure 5: Greater values of 9 co i result in dust echoes that are brighter at intermediate times 
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because more light is available to be scattered by dust. 

The smaller values of R that we consider in this paper permit two simplifications: (1) We 
have assumed that dust within Oj e t is destroyed to well beyond R (which is equivalent to setting 
/ = 1 in Reichart 2001); and (2) We can limit the discussion to light scattered from the optical 
flash, which we take to be light from either internal shocks (e.g., Vestrand et al. 2005) or the 
reverse shock of the afterglow: Reichart (2001) showed that light scattered from the forward shock 
of the afterglow can only outshine the afterglow for large values of R. 

Consider the case of an optical flash of the form: 

where F Vt $ is the brightness of the optical flash and to is the timescale of the optical flash, after 
which it fades as t~ 2 . The spectral fiuence of the optical flash is then 2F u fltQ and by Equation 24 
of Reichart (2001) the peak brightness of the dust echo is given by: 

F »-™~ 0.0045 (Z^V_« rV <5J), (5) 



F v>0 ' V 1 + z J \3 x 10 14 cm J V 30 sec 

where F Ut oE is the peak brightness of the dust echo, t v !\ +z -\ is the optical depth through the dust 
distribution at frequency v(l + z), and z is the redshift. 2 By Equation 5 of Reichart (2001), the 
peak time of the dust echo is given by: 

where tpE is the peak time of the dust echo. For these parameter values, 
\og(F v ^DE / F v fi) / \og{tDE /to) ~ —1-9 and consequently the dust echo would outshine the 
optical flash beginning at toE ~ 8 min, which would result in a «0.9 mag bump in the light 
curve at this time, if the forward shock is not yet a major contaminant. Furthermore, the dust 
echo would be even more prominent at bluer wavelengths, since t v i 1+z \ is an increasing function 
of v. Greater values of T u ti +Z ^, R, z, and especially 6j e t, and lesser values of to, might also make 
the dust echo more prominent, in part by delaying it. However, the forward shock is more of a 
contaminant at later times as well. Even so, for T u n +Z \, R, and 8j et three times larger and to 
three times smaller, \og(F v ^E / F v $) / \og(toE /to) ~ —0.9, which means that the dust echo would 
outshine the forward shock beginning at tu-E ~ 11 hr as long as the optical flash outshines the 
forward shock at t ~ 10 sec (for a forward shock that fades faster than i -0 - 9 ). 

However, these estimates for the scattering of light by grains of dust must also be taken with 
a few grains of salt. The model of Reichart (2001) makes a number of limiting assumptions: (1) 



2 Equation 24 of Reichart (2001) has a factor of e r "( 1 +*) that we have dropped here, because Equation 24 is an 
approximation that does not hold around t « toE when t„( 1+z ) > 1. However, from Figure 7 of Reichart (2001) it is 
clear that dropping this factor approximates the correct effect. 
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The viewing angle is assumed to be zero; (2) Light that scatters is assumed to not scatter again 
into the line of sight; and (3) The radius at which the optical flash occurs is assumed to be much 
less than R. 

The effect of a non-zero viewing angle will be to stretch out the turn-on time of the dust 
echo, from t ~ Q.Ude to a delayed peak time of t ~ 2.7toE for the typical viewing angle (Reichart 
2001). The light curve and color history of this transition is difficult to predict without performing 
numerical simulations. The scattering assumption should be fine for Tj,( 1+2 ) < 1 but will become 
increasingly less accurate at greater optical depths. 

Although internal-shock optical flashes should occur at smaller radii and not be a problem, 
reverse-shock optical flashes should begin at the deceleration radius, which for a wind-swept 
environment is given by: 

r]£ v. 



where r\ is a constant of proportionality that is equal to 9/2 in the case of the Blandford-McKee 
solution, £ is the isotropic-equivalent energy release of the explosion, v w is the velocity of the 
progenitor's wind, M is the mass-loss rate of the progenitor, and Tq is the initial Lorentz factor 
of the ejecta (e.g., Kumar k Panaitescu 2003). For 77 = 9/2, £ = 10 52 erg, v w = 2,000 km/s, 
M = 5 x 10~ 5 M /yr, and T = 300, Td = 3 x 10 cm. Since this is considerably less than the 
values of R that Reichart (2001) considered, he was able to safely ignore this offset from the center 
of the dust distribution in his model, which assumes spherical symmetry (except along the jet axis 
where dust is destroyed). However, now this case is marginal: For larger values of £ in particular, 
this model will no longer apply and light curves and color histories will be difficult to predict 
without performing numerical simulations. However, given the standard-energy result (e.g., Frail 
et al. 2001; Bloom et al. 2003) larger values of £ imply smaller values of 9j et , the dust echoes for 
which would be difficult to detect anyway (Equation 6). 

However, the model of Reichart (2001) might still apply to reverse-shock optical flashes, even 
for large values of £. For example, if the central engine is long lived only a fraction of this energy 
might be released initially, in which case the optical flash would occur at a smaller distance. 
Also, R might be greater than what we expect from late-type WC stars: If the progenitor's final 
years are accompanied by greater ultraviolet output, dust will be sublimated to greater distances. 
Furthermore, greater charging of carbon atoms might prevent grain formation altogether (e.g., 
Williams, van der Hucht & The 1987), in which case R would be set by the wind velocity and the 
time since last significant grain formation, possibly increasing R by a factor of ten or more. 



4. Conclusions and PROMPT 



We have revisited the possibility of GRB dust echoes in light of expectations for a WC star 
progenitor. We find that if the optical flash is sufficiently bright and the jet sufficiently wide, 
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GRB afterglows may be accompanied by detectable chromatic variations, not on a timescale of 
weeks after the burst but on a timescale of minutes to hours. From these signatures, such model 
parameters as the inner radius of the dust distribution, the initial opening angle of the jet, etc., 
may be deduced. 

However, in order for early-time dust echoes to be identified and studied, three conditions 
must first be satisfied: (1) GRBs must be localized rapidly, preferably within tens of seconds; (2) 
Optical/near- infrared telescopes must respond on a similar timescale; and (3) These telescopes 
must observe in many wavelength bands and preferably simultaneously, else it may be impossible 
to disentangle SFDs from temporal variability, particularly at early times. Condition (1) is being 
satisfied by HETE-2, Integral, and now Swift, about once every fourth day. Condition (2) is being 
satisfied by a wide variety of robotic telescopes. 

The Panchromatic Robotic Optical Monitoring and Polarimetry Telescopes (PROMPT) that 
the University of North Carolina at Chapel Hill is building at the Cerro Tololo Inter- American 
Observatory will additionally satisfy condition (3). When completed in late-2005, PROMPT will 
consist of six 0.41-meter Ritchey-Chrtien telescopes on rapidly slewing (9° /sec) Paramount ME 
mounts and will image GRB afterglows in nine wavelength bands - u'g'r'Ri'z'YJH - six of them 
simultaneously. (The R telescope will additionally measure the polarization of the afterglow.) In 
addition to measuring redshifts by dropout and extinction curves in great detail, especially when 
combined with Swift UVOT and XRT measurements, PROMPT'S rapid slewing and simultaneous 
multi-wavelength imaging design will make it ideal for identifying and studying dust echoes at 
early times. 

DER very gratefully acknowledges support from NSF's MRI, CAREER, PREST, and REU 
programs, NASA's APRA, Swift GI, and IDEAS programs, Dudley Observatory's Ernest F. 
Fullam Award, and especially Leonard Goodman and Henry Cox. 
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Fig. 1. — Top panel: Dust and light geometries and model parameters. Dust does not form within 
radius R of the progenitor and is destroyed within half-angle 9j e t of the jet axis by bipolar jets 
(e.g., Waxman & Draine 2000; Fruchter, Krolik & Rhoads 2001). Exterior to the solid curve, the 
dust density n(r) oc r -2 . The optical depth through this dust distribution is r v = n(R)a u R, where 
a v is the total (absorption plus scattering) cross section of the dust grains at frequency v. A pulse 
of light is assumed to be emitted within half-angle 9 co i of the jet axis (Reichart 2001). Bottom 
panel: Curves of constant arrival time t = 0.01, 0.03, 0.1, 0.3, 1 and 3 in units of R(l + z)/c. Solid 
portions mark illuminated dust. 

Fig. 2. — Light curves (top panel) and color histories (bottom panel) of dust echoes for R(l + z) = 
10 14 , 3 x 10 14 ' 10 15 , 3 x 10 15 , and 10 16 cm, 9 jet = 10°, t u(1+z) = 1, and 9 col = 1.19 jet . The light 
curves are normalized to the fluence of the optical flash divided by one day. The parameter a V) DE{t) 
is the spectral index of the dust echo, a u is the spectral index of the optical flash, and oc v / 1+z -\ is 
given by t u ^ +z ^ cc \u{1 + z)] a "( 1+z '> . The parameter a u is likely 1/3, —1/2, or — (p — l)/2, where p 
is the power- law index of the electron-energy distribution (e.g., Sari, Piran & Narayan 1998), and 
1 £ a„(i+*) £ 1.6 (Reichart 2001). 

Fig. 3. — Light curves (top panel) and color histories (bottom panel) of dust echoes for 9j e t = 2.5°, 
5°, 10°, 20°, and 40°, R(l + z) = 10 15 cm, t u{1+z) = 1, and 9 col = 1.19 jet . 

Fig. 4. — Light curves (top panel) and color histories (bottom panel) of dust echoes for Tj,(i +2 ) = 0.1, 
0.3, 1, 3, and 10, and R(l + z) = 10 15 cm, 9 jet = 10°, and 9 coi = 1.19 jet . 

Fig. 5. — Light curves (top panel) and color histories (bottom panel) of dust echoes for 9 co i/9j e t = 
1.01, 1.03, 1.1, 1.4, and 2.7, R(l + z) = 10 15 cm, 9 jet = 10°, and t u(1+z) = 1. 
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